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INTRODUCTION 

. The  objective  of  the  MM&T  Project  6550  was  to  perform  engineering  studies  and 

surveys  of  industry  in  order  to  assure  that  the  methods  of  shell  production  in  modernized 
Army  plants  were  consistent  with  the  latest  state  of  art  in  manufacturing  technology. 

This  task  dealt  with  the  modeling  and  computer  simulation  of  the  drawing  operation  in 
shell  manufacture.  Specifically,  for  shell  drawing,  the  task  objectives  were:  (a)  to 
provide  analytical  models  describing  the  forming  of  artillery  shell  and  cartridge  cases, 
and  (b)  to  optimize  the  modeled  process.  The  work  performed  under  this  task  is  pre- 
j sented  in  this  report. 

In  current  shell  manufacturing  practice,  the  shell  is  formed  to  have  an  almost 
finished  internal  cavity  surface  with  the  finish  machining  largely  relegated  to  the  out- 
side surface  of  the  piece.  Nosing  then  occurs  after  rough  machining  is  completed. 

There  are  three  basic  methods  of  introducing  a shell  cavity  into  the  slug  or  billet:  (1) 

Hot  piercing  followed  by  hot  drawing,  (2)  Cold  extrusion  and  drawing  and  (3)  Hot  cup- 
ping followed  by  cold  drawing.  In  each  of  these  methods  of  initial  shell  forming,  the 
drawing  process,  whether  hot  or  cold,  is  the  essential  basic  step  of  the  several  form- 
ing operations  involved.  Consequently,  due  to  its  basic  technological  importance,  draw- 
ing was  selected  as  the  first  operation  to  be  modeled. 

This  study  was  also  undertaken  because  of  a lack  of  data  in  manufacturing  using 
mechanical  presses.  Projected  shell  lines  will  use  mechanical  presses,  or  a combina- 
tion of  mechanical  and  hydraulic  presses,  rather  than  the  conventional  purely  hydraulic 
systems.  Data  on  draw  feasibility  and  process  design  with  a mechanical  press  for 
artillery  shell  is  very  limited.  This  sparsity  of  data  on  processing  with  mechanical 
presses  at  elevated  speeds  presents  design,  planning  and  commonality  problems  for 
both  current  and  projected  production  lines.  The  models  of  the  drawing  process,  as 
documented  in  this  report,  can  provide  reliable  guidelines  for  process  design  with 
mechanical  press  systems.  This  study  on  shell  drawing  supports  all  current  and  pro- 
jected 105mm  and  155mm  lines.  The  computer  programs  generated  under  this  task 
are,  however,  general  enough  to  support  any  size  shell. 

The  objectives  of  this  task  were  met  by  a combination  of  in-house  and  contractural 
effort  with  Battelle-Columbus  Laboratories.  Battelle1  developed  the  math  models 
associated  with  the  curved  die  profiles,  optimization  and  a portion  of  the  heat  transfer 
work.  These  efforts  have  been  thoroughly  documented  in  a report1  which  shows  the 
details  of  the  models  and  the  computer  programs.  Part  of  the  Frankford  effort  con- 
sisted of  developing  math  models  associated  with  failure  during  drawing,  heat  transfer 
and  conventional  tapered  die  profiles.  The  present  report  unifies  these  two  develop- 
ment efforts  by  presenting  results  drawn  from  both  sources.  All  models  generated 
, under  this  task  have  been  made  functional  on  the  CDC  6000  series  computers  at 


(1)_  . 

References  are  given  at  the  end  of  this  report. 
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Picatinny  Arsenal  via  remote  entry  from  Frankford  Arsenal.  All  the  numerical  results 
presented  in  this  report  were  generated  at  Frankford  Arsenal. 

As  a result  of  this  study  a new  curved  die  configuration  was  designed  which  requires 
a minimum  energy  expenditure  from  the  press  for  any  specified  reduction,  either  hot 
or  cold.  The  new  die  should  require  10-15%  less  energy  than  a conventional  straight 
profile  ring  performing  the  same  reduction.  The  effects  of  processing  conditions  such 
as  temperature,  punch  speed,  friction  and  reduction  for  the  new  configuration  were 
studied  in  order  to  present  criteria  which  avoid  punch  through  failure.  Concrete  ex- 
ample designs  were  performed  for  the  cold  drawing  of  a 105mm  shell  forging  using  a 
variety  of  current  steels.  It  was  found  that  40-42%  reduction  in  a single  pass,  cold 
operation  is  the  upper  limit  of  draw  feasibility  for  a 105mm  shell.  Treated  also  in  this 
report  were  conventional  tapered  dies.  Design  curves  covering  a wide  range  of  pro  - 
cessing  variables  including  failure  were  prepared.  In  general,  for  boih  curved  and 
tapered  dies,  it  was  established  that  the  temperature  in  hot  forming  plays  a key  role 
in  limiting  potential  operating  press  speed. 

The  remainder  of  this  report  is  organized  as  follows:  specific  data  dealing  with 
optimal  curved  dies  and  conventional  tapered  dies  are  presented  in  the  Simulation  Re- 
sults section.  Interpretation  of  the  data  and  conclusions  are  presentee  in  the  Summary 
and  Discussion.  The  mathematical  models  used  have  been  placed  sep; irately  in  Appendix 
A and  B,  rather  than  in  the  text  proper.  It  is  not  necessary  for  users  of  the  data  in 
this  report  to  refer  to  these  Appendicies.  Previous  models,  a literati  re  survey  and 
assessments  of  the  state  of  the  art  in  simulation  have,  because  of  thei;-  length,  been 
placed  in  Appendix  A. 
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I OPTIMAL  CURVED  DIES 
(1) 

On  the  basis  of  the  Battelle  study,  it  was  found  that  curved  die  profiles  could  be 
designed  which  were  optimal.  Optimal  in  this  context  has  the  following  meaning:  for 
a given  reduction,  material  and  lubrication  (friction),  condition  there  is  a curved  die 
profile  which,  when  compared  to  a conventional,  tapered  profile  is  superior  in  that: 

a.  Less  energy  is  required  to  form, 

■ b.  The  flow  of  metal  is  such  that  the  strains  at  the  exit  are  uniform  through  the 

shell  wall. 

The  first  factor  implies  a potential  cost  reduction  by  conserving  press  power.  The 
second  factor  implies  a reduction  in  liklihood  of  material  failure  or  "punch  through"  by 
the  mandrel.  Experimental  studies  indicate  that  the  product,  due  to  the  radial  homo- 
geneity of  the  strains  developed,  will  also  show  improved  metallurgical  properties. 

This  report  extends  the  above  findings  by  considering  the  influence  of  various 
processing  variables  and  materials  on  the  design  of  optimal,  curved  dies.  The  com- 
puter program  CDVEL  was  used.  An  extended  discussion  of  the  mathematical  model 
in  CDVEL  is  to  be  found  in  Reference  1,  and  is  also  summarized  in  section  II.  A of 
Appendix  A of  this  report. 

Figure  1 shows  the  problem  under  consideration.  A shell  is  being  drawn  through 
a die  (cross  hatched)  whose  inner  surface  is  doubly  curved  by  two  circular  arcs  of 
radii  R1  and  R„.  A punch  of  radius  rQ  is  moving  the  material  from  left  to  right  with 
an  entry  speed  v0  and  an  exit  speed  of  Ve.  The  basic  problem  is  to  determine  the 
profile  (i.  e. , Rj  and  R2)  in  the  optimal  manner  described  above. 

k 

The  results  to  be  presented  are,  as  in  Reference  1,  for  the  40%  reduction  of  a 
105mm  shell  in  a single  pass  with  an  untapered  punch.  The  product  thus  has  a constant 
wall  thickness.  The  baseline  conditions  for  cold  and  hot  drawing  used  here  are  in 
Table  1.  Departures  from  these  values  are  indicated  in  the  text  as  necessary. 

The  effect  of  using  various  commonly  used  shell  steels  on  optimal  shape  in  a cold 
40%  draw  are  shown  in  Figure  2.  This  graph  represents  the  actual  contour  (in  mm)  of 
the  optimal  curved  die  for  each  steel.  The  radial  coordinate  (R(z)  of  Figure  1)  is 
plotted  as  a function  of  the  axial  coordinate  (z  of  Figure  1).  The  actual  data  points, 
which  are  automatically  printed  by  the  computer,  are  given  in  Table  2.  (Except  for 
two  steels,  52100  and  1095,  the  material  properties  are  taken  at  20° C - cold.  The 
52100  and  1095  runs  were  made  at  750°  C,  the  lowest  temperature  for  which  material 
■ properties  were  available.  As  can  be  immediately  seen,  the  optimal  die  profile  is 

essentially  the  same  for  each  steel  except  for  the  4337).  Thus,  the  important  con- 
clusion: For  cold  and  moderately  warm  drawing  at  moderate  speeds  the  optimal  die 
profile  is  essentially  the  same  for  all  steels. 


Table  1.  Conditions  Used  for  Theoretical  Determination  Of  The  Die 
Shape  And  For  The  Prediction  Of  Process  Variables 


Billet  Material  = 

Billet  Outside  Diameter 
Billet  Inside  Diameter 

Punch  Diameter  = 

Reduction  in  Area  of  Cross  Section 
Punch  Speed  = 

Friction  Shear  Factor  mj  (die/shell) 

Friction  Shear  Factor  m2  (die/mandrel)  = 

Billet  Temperature  = 


AISI  1045  Steel(*) 
140  mm 
86  mm 
86  mm 
40  Percent 
100  mm/sec  (cold) 
1000  mm/sec  (hot) 
0.1  (cold) 

0.25  (hot) 

0.2  (cold) 

0.35  (hot) 

20°  C (cold) 

1000°  C (hot) 


(*)The  flow  stress  data  for  this  material  were  obtained  from  Reference  29. 


Table  3 shows  the  drawing  data  for  the  above  steels.  Examining  the  specific  draw 
stress,  which  represents  failure  by  punch  through  if  it  exceeds  unity  (Appendix  A): 

(the  1095  and  52100  steels  are  totally  inadequate  at  750°  C;  with  40%  reduction  failure 
will  certainly  occur.  The  same  steels  in  the  cold  condition  may  be  adequate  at  a 40% 
reduction. ) The  4337  steel  is  marginal  for  failure;  the  remaining  steels  should  not 
experience  failure. 

The  energy  rates  for  the  failed  steels  are  much  lower  than  for  the  safe  materials. 
Similar  commentary  applies  to  the  drawing  load  and  the  punch  load.  These  facts  are 
explained  by  the  fact  that  the  yield  strength  has  been  degraded  in  the  two  metals  which 
failed.  Thus  a much  lower  energy  and  load  system  resulted.  Thus  the  definition: 

In  designing  an  optimal  die  for  a given  material,  temperature  mandrel  speed  and 
friction  the  safe  operating  condition  is  such  that  the  specific  draw  stress  (line  stress 
ratio)  must  be  less  than  1. 

In  the  Battelle^  report,  it  was  shown  that,  in  general,  the  influence  of  work  or 
strain  hardening,  temperature  and  mandrel  speed  does  not  appreciably  effect  the  shape 
of  the  optimal  die.  However,  as  noted  above,  failure  is  heavily  dependent  on  these 
factors.  Shown  in  Figures  3 and  4 are  the  drawing  limit  diagrams  for  1045  steel  for  a 
wide  range  of  frictional  conditions.  The  ordinate  is  the  specific  draw  stress,  the 
absissa  reduction  in  area  in  a single  pass.  In  Figure  3 the  mandrel  friction  was  con- 
stant and  the  shell  die  friction  varied.  The  reverse  is  true  of  Figure  4.  It  may  be 
noted  that:  Regardless  of  frictional  conditions,  40-42%  is  the  largest  possible  single 
pass  cold  draw  reduction  on  a 105mm  shell'  Also:  Reducing  the  friction  between  shell 
and  die  increases  the  maximum  reduction  which  can  be  taken.  However,  increasing  the 
friction  between  shell  and  mandrel  increases  the  maximum  reduction  which  can  be  taken. 
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Table  2.  R And  Z Coordinates  For  Optimal  Die  Contours  Shown  in  Figure  2. 
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REDUCTION  IN  AREA,  PERCENT 

Figure  3.  Drawing  Limit  Diagram  for  a Constant  Friction  at  the  Punch-Material 
Interface  (RQ  = 70  mm,  r0  = 43  mm.  Temperature  = 20° C,  Curved- 
Straight-Curved  Die  Angle  = 17.  8 Degrees) 


1.2  mi  = 0.10 
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REDUCTION  IN  AREA,  PERCENT 

Figure  4.  Drawing  Limit  Diagram  for  a Constant  Friction  at  the  Die-Material 
Interface  (R0  = 70  mm,  rQ  = 43  mm,  Temperature  = 20°  C,  Curved- 
Straight-Curved  Die  Angle  = 17.  8 Degrees) 
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The  effect  of  friction  on  shape  is  shown  in  Figure  5.  Shown  are  the  optimal  40% 
reduction  contours  for  1045  steel  under  extremes  of  frictional  conditions.  Curve  1 
has  the  base  line  friction  conditions  of  Figure  2.  Curve  2 represents  the  desirable 
quality  of  decreasing  the  shell  die  friction.  Curve  3 represents  the  desirable  quality 
of  increasing  the  friction.  The  conclusion  was  that:  As  in  the  case  of  strain-harden- 
ing and  temperature,  friction  does  not  significantly  effect  the  optimal  die  contour. 

Figures  6-14  show  the  effect  of  reduction  in  area  on  die  contours  for  the  drawing 
of  a 1045  steel  105mm  billet.  Each  plot  represents  the  optimal  die  contour.  The 
computed  coordinates  are  in  Table  4,  and  the  failure  design  data  is  in  Table  5.  The 
specific  draw  stress  data  of  Table  5 shows,  as  mentioned  above,  that  40-42%  is  the 
maximum  possible  reduction  in  a cold  draw.  It  is  also  noted  that  as  the  areal  reduction 
is  increased,  the  die  length,  drawing  load  and  punch  loads  also  increase.  The  drawing 
load  is  the  total  load  which  must  be  carried  in  the  shell  wall,  and  the  punch  load  is  dis- 
tributed over  the  punch  face.  The  effect  of  increased  punch  load  will,  of  course,  be  to 
reduce  tool  life.  The  strain  differential  column  of  Table  5 represents  the  differential 
in  plastic  strain  between  the  outermost  and  innermost  fibers  in  the  shell  wall  at  the 
exit.  It  is  therefore  a measure  of  the  strain  uniformity  in  the  product  wall.  As  can 
be  noted,  an  increased  areal  reduction  leads  to  greater  lack  of  product  homogeneity. 
However,  as  shown  in  the  Richmond/Devenpeck  study,  for  any  reduction  the  stream- 
lining effect  of  curved  dies  will  lead  to  a more  uniform  product. 
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Figure 


Figure  6.  Optimal  Die  Contours  for  1045  Steel  in  a Single  Pass  Draw  (See  Table  1)  with  Various  Reductions 


Figure  8.  Optimal 


Figure  10.  Optimal  Die  Contours  for  1045  Steel  in  a Single  Pass  Draw  (See  Table  1)  with  Various  Reductions. 


Figure  11.  Optimal  Die  Contours  for  1045  Steel  in  a Single  Pass  Draw  (See  Table  1)  with  Various  Reductions. 


Figure  12.  Optimal  Die  Contours  for  1045  Steel  in  a Single  Pass  (See  Table  1)  with  Various  Reductions 


Table  4.  R and  Z Coordinates  for  Optimal  Die  Contours  Shown  in  Figures  6-14 

OPTIMAL  DIE  CONTOURS 


29.  3% 

32.  0% 

36. 

0% 

40.  0% 

Z 

r 

Z 

r 

Z 

r 

Z 

r 

0 

70.  00 

0 

70.  000 

0 

70. 000 

0 

70.  000 

5.  946 

69.  732 

6.  22 

69.  707 

6.608 

69.  669 

6.968 

69. 625 

10.901 

69.  096 

11.403 

69.  01 

12. 115 

68.  882 

12.  775 

68.  732 

14.866 

68.  309 

15. 549 

68. 148 

16.521 

67. 905 

17.420 

67.622 

17. 839 

67. 551 

18.  659 

67. 316 

19.  825 

66. 962 

20.905 

66.  546 

19. 821 

66.  963 

20.732 

66.  669 

22.  027 

66. 227 

23.  227 

65.708 

20. 812 

66.  643 

21. 769 

66.318 

23. 129 

65. 828 

24.389 

65. 251 

21.802 

66. 324 

22. 804 

65. 967 

24. 229 

65.429 

25.511 

64. 809 

23. 782 

65. 736 

24. 876 

65. 301 

26.  429 

64. 695 

27.  756 

63.  998 

26.752 

64. 979 

27.982 

64. 490 

29.730 

63. 753 

31. 124 

62.  959 

30.712 

64. 183 

32. 124 

63.628 

34. 131 

62.777 

35.614 

61. 886 

35.661 

63. 557 

37. 302 

62.  933 

39.632 

61.990 

41.227 

61.022 

41.601 

63. 290 

43.515 

62. 640 

46.233 

61.660 

47.963 

60.660 

41. 

0% 

42.  0% 

44.  0% 

48.  0% 

Z 

r 

Z 

r 

Z 

r 

Z 

r 

0 

70.  000 

0 

70. 000 

0 

70.  000 

0 

70. 000 

7.034 

69.616 

7.  066 

69.  610 

7.309 

69.584 

7.359 

69.  549 

12. 896 

68.699 

12.954 

68.  680 

13.400 

68.591 

13.492 

68.471 

17.586 

67.558 

17.664 

67.  523 

18.  273 

67. 355 

18.  398 

67. 124 

21. 103 

66.453 

21. 197 

66.401 

21.  928 

66.  154 

22. 077 

65. 813 

23.  448 

65.591 

23.  552 

65.  525 

24.  364 

65.216 

24.  530 

64. 784 

24.620 

65. 120 

24.729 

65.  047 

25.583 

64.704 

25.757 

64.  221 

25.752 

64. 606 

25.  894 

64.575 

26.745 

64. 215 

26.948 

63.  675 

28. 015 

63. 834 

28.  223 

63.709 

29.  071 

63.  320 

29.33 

62.  676 

31.410 

62. 767 

31.716 

62. 600 

32.559 

62. 175 

32.903 

61.402 

35. 937 

61.666 

36. 374 

61.455 

37. 210 

60. 994 

37.667 

60. 095 

11.596 

60.781 

42.  197 

60.  536 

43.  023 

60.  047 

43.612 

58.  048 

18.387 

60.410 

49. 184 

60. 150 

49.  999 

59.650 

50.768 

58.610 

51. 

Z 

1% 

r 

51.1  (Continued) 
Z r 

0 

70. 000 

29.  804 

62.  144 

7.480 

69.519 

33.428 

60.  779 

13.  713 

68.  367 

38.  260 

59.  382 

18.700 

66.926 

44.300 

58.  266 

22.  440 

65.517 

51.548 

57. 800 

24. 933 

64.411 

26. 179 

63,  804 

27. 387 

63. 216 
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n CONVENTIONAL  TAPERED  DIES 


In  order  to  address  the  problems  of  simulation  and  failure  in  conventional  dies, 
a computer  program  SHELL  has  been  developed.  A mathematical  summary  is  given 
in  Appendix  A.  The  purpose  of  this  section  is  to  present  guidelines  for  designing  tap- 
ered dies  which  will  not  permit  material  failure  (punch-through). 

Table  6 indicates  the  baseline  process  variables  being  used.  Changes  in  these 
quantities  are  noted  when  appropriate.  Figures  15  and  16  are  drawing  limit  diagrams 
for  a cold,  40%  reduction  of  AISI  1045  billet  with  various  die  angles.  These  figures 
plot  specific  draw  stress  vs  die  semi  angle.  As  in  curved  dies,  when  the  specific  draw 
stress  is  unity  failure  occurs.  Figure  15  shows  that,  for  any  level  of  friction  between 
shell  and  die,  increasing  the  die  angle  avoids  failure.  Decreasing  the  shell-die  fric- 
tion enhances  the  process.  Figure  15  also  demonstrates  that  for  any  friction  level 
there  is  a minimum  angle  to  avoid  failure.  Figure  16  shows  the  effect  of  varying  ^2, 
the  shell  mandrel  friction.  Increasing  this  friction  enhances  the  process  in  that  fail- 
ure cannot  occur.  Noting  that  drawing  with  zero  or  negative  draw  stress  that  there  is 
a minimum  angle  below  which  drawing  cannot  occur. 


Table  6.  Conditions  Used  for  Theoretical  Determination  of  Process 
Variables  in  Conventional  Dies 


Billet  Material 

= 

AISI  1045  Steel  (*) 

Billet  Outside  Diameter 

= 

140  mm  (**) 

Billet  Inside  Diameter 

= 

86  mm  (**) 

Punch  Diameter 

= 

86  mm  (**) 

Reduction  in  Area  of  Cross  Section 

= 

40  Percent  (**) 

Punch  Speed 

= 

100  mm/sec  [cold] 

= 

1000  mm/sec  [hot] 

Friction  Coefficient,  (shell/mandrel) 

= 

0. 1 [cold] 

= 

0.25  [hot] 

Friction  Coefficient,  /j  2 (shell/die) 

= 

0.2  [cold] 

= 

0.  35  [hot] 

Billet  Temperature 

= 

20°  C [cold] 

= 

1000°  C [hot] 

Die  Half  Angle  a 

= 

Variable 

(*)  The  flow  stress  data  for  this  material  were  obtained  from 
Reference  29. 

(**)These  vary  with  reduction  in  area. 
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Processing  Conditi 


Figure  17  and  18  present  an  analagous  study  to  the  one  above,  except  that  drawing 
was  performed  at  1000°  C,  Comparison  of  Figures  15  and  17  shows  that  hot  drawing  is 
a similar  process  except  that  the  sp.  draw  stress  is  lower.  A comparison  of  Figures 
16  and  18  shows  that  in  both  hot  and  cold  drawing,  increased  mandrel  friction  is  bene- 
ficial to  avoiding  failure. 

The  effect  of  the  punch  velocity  on  cold  drawing  is  shown  in  Figure  19.  The  plot 
of  specific  draw  stress  vs  semiangle  at  six  punch  velocities  is  shown.  All  curves 
coincide.  This  may  be  explained  in  terms  of  the  materials  properties  in  Table  7. 

Table  7 shows  the  dependence  of  the  metals  used  in  this  on  strain  hardening  and  strain 
rate  sensitivity.  Strain  hardening  is  covered  by  assuming  that  the  equivalent  plastic 
stress  (o)  is  proportional  (k)  to  the  plastic  strain  raised  to  a power  (n).  This  is  seen 
to  be  true  for  cold  to  warm  (400-500°  C)  temperatures.  In  this  range  the  value  of  k 
and  n is  somewhat  influenced  by  strain  rate.  For  temperatures  above  600°  C,  the  above 
law  is  not  valid.  In  this  range  a better  assumption  is  that  the  stress  (o)  is  proporational 
(c)  to  the  strain  rate  (f).  The  values  of  c and  n depend  upon  the  level  of  strain.  Since 
the  velocity  of  the  punch  defines  the  level  of  strain  rate  it  should  be  anticipated  that, 
at  higher  temperatures,  this  factor  is  important,  and,  conversely  at  lower  tempera- 
tures unimportant. 

Figure  19  indicates  the  general  conclusion  that  the  punch  velocity  does  not  in- 
fluence cold  drawing.  ' " ' 


Figure  20  shows  the  influence  of  punch  velocity  in  a 40%  hot  (1000°  C)  draw  with 
moderate  friction  and  punch  velocities  from  100  to  1000  mm/s.  The  trend  is  the  same 
as  in  cold  drawing. 

An  additional  effect  of  temperature  for  a low  speed  draw  is  illustrated  in  Figure 
21.  Shown  is  the  drawing  limit  diagram  for  a 100  mm/s  draw  over  a range  of  temper- 
atures from  50°  C to  1000°  C.  It  is  seen  that  up  to  about  500°  C,  the  cold-warm  range, 
the  temperature  has  no  effect,  and  the  specific  draw  stress  is  constant  with  increasing 
half  angle.  Between  500°  C and  1000°  C the  specific  draw  stress  is  lowered  by  about 
10%,  and  the  drawing  curve  is  parallel  to  the  cold  curve.  Beyond  1000°  (not  shown) 
the  behavior  is  nonlinear.  As  the  molten  state  is  approached  the  draw  stress  drops 
and,  as  should  be  expected,  drawing  is  not  possible. 

A corresponding  study  for  a high  punch  speed  is  shown  in  Figure  22.  The  punch 
speed  is  1000  mm/s,  ten  times  that  of  Figure  21.  Up  to  about  1000°  C the  material 
response  is  essentially  the  same  as  in  drawing  at  100  mm/sec. 

A final  comparison  between  hot  and  cold  drawing  is  displayed  in  Figures  23  and 
24.  As  is  shown  in  Figure  23  when  the  coefficients  of  friction  are  equal,  irrespective 
of  their  value,  the  curve  is  a horizontal  line.  This  is  the  one  circumstance  in  cold 
drawing  under  which  the  die  angle  is  not  important.  Figure  24  shows  a similar  study 
at  1000°  C as  can  be  seen  as  the  same  as  for  a cold  draw. 
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FAILURE 


Specific  Draw  Stress  vs.  Die  Half  Angle  for  Various  Processing  Conditions 


J.8.  Specific  Draw  Stre 


Specific  Dr^nv  Stress  vs.  Die  Half  Angle  for  Various  Processing  Conditions 


Table  7.  Material  Properties 


I.  Cold  (20°  C)  Drawing 
o = K (e)n 


Steel 

Strain  Rate,  l/s 

Kd  103  psi 

n 

1015 

1.6 

113.8 

.10 

1045 

1.6 

147.9 

.11 

1045 

1.5 

137.9 

.14 

II.  Hot  ( 600°  C)  Drawing 
a = C (ihm 

Steel  Strain!  600°  800°  1000°  1200°-  Temp°C 


c 

m 

c 

m 

c 

m 

c 

m 

1015 

.2 

36.8 

.112 

1015 

.25 

19.9 

. 105 

17.0 

.045 

7.2 

.137 

1015 

.4 

40.  6 

.131 

1015 

.5 

21.5 

.104 

18.8 

.058 

6.8 

.169 

1015 

.6 

40.0 

.121 

1015 

.7 

39.5 

.114 

21.1 

.109 

18.3 

.068 

5.7 

.181 

900°  1000°  1100°  1200°-  Temp"  C 


c 

m 

c 

m 

c 

m 

c 

m 

1045 

.05 

25.4 

.080 

15.1 

.089 

11.2 

.100 

8.0 

.175 

1045 

.10 

28.9 

.082 

18.8 

.103 

13.5 

.125 

9.4 

. 168 

1045 

.20 

33.3 

.086 

22.8 

.108 

15.4 

.128 

10.5 

.167 

1045 

.30 

35.4 

.083 

24.6 

. 110 

15.8 

.162 

10.8 

.180 

1045 

.40 

35.4 

.105 

24.7 

.134 

15.5 

. 173 

10.8 

.188 

1095 

.10 

18.3 

. 146 

13.9 

.143 

9.8 

.159 

7.1 

.184 

1095 

.30 

21.9 

.133 

16.6 

.132 

11.7 

.147 

8.0 

.183 

1095 

.50 

21.8 

.130 

15.7 

.151 

10.  6 

.176 

7.3 

.209 

1095 

.70 

21.0 

.128 

13.6 

. 179 

9.7 

.191 

6.5 

.232 

4337 

.10 

22. 1 

.080 

16.6 

.109 

12.1 

.115 

8.2 

. 165 

4337 

.30 

28.1 

.077 

20.8 

.098 

15.0 

.111 

10.7 

.138 

4337 

.50 

29.2 

.075 

21.8 

.096 

15.7 

.112 

11.3 

.133 

4337 

.70 

28.1 

.080 

21.3 

. 102 

15.5 

.122 

11.3 

.135 

52100 

.10 

20.9 

.123 

14.3 

.146 

9.5 

.169 

6.7 

.203 

52100 

.30 

25.5 

.107 

17.7 

.127 

12.0 

.143 

8.3 

.171 

52100 

.50 

25.9 

.107 

17.7 

.129 

12.3 

.143 

8.3 

.17,8 

52100 

.70 

23.3 

.131 

16.8 

.134 

12.0 

.148 

7.7 

.192 
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Figure  23.  Specific  Draw  Stress  vs.  Die  Half  Angle  for  Various  Processing  Conditions. 


It  should  be  noted  with  caution  that  much  of  the  above  study  was  performed  with 
the  two  coefficients  of  friction  equal.  In  practice  this  equality  condition  is  seldom 
seen.  The  striking  similarity  shown  between  cold  and  hot  forming  will  hold  only  in 
this  restricted  condition.  As  can  be  noted  from  Figure  16  and  18,  when  the  two  co- 
efficients differ  this  results  in  a larger  difference  between  hot  and  cold  drawing. 

This  difference  manifests  itself  as  a reduced  stress  in  hot  drawing. 

Ill  HEAT  TRANSFER  IN  DRAWING 

As  mentioned  in  the  Battelle  study  and  in  the  Appendix  to  this  report,  a study  of 
the  effects  of  heat  conduction  during  hot  drawing  was  undertaken.  In  drawing  both 
plastic  flow  and  friction  contribute  to  heat  generation.  About  90-95%  of  the  mechanical 
press  energy  expended  in  forming  will  appear  as  heat.  As  was  seen  in  Section  n above, 
temperature  at  higher  speeds  is  the  principal  factor  in  failure.  The  mathematical 
model  used  to  predict  the  conduction  of  heat  is  described  in  the  Battelle report.  The 
computer  programs  resulting  from  that  report  were  used  to  achieve  the  results  given 
below.  The  conditions  used  were  for  40%  hot  reduction  of  1045  steel  as  shown  in 
Table  1. 

Figure  25  shows  the  isotherms  (lines  of  constant  temperature)  when  steady  state 
conditions  are  reached.  The  initial  billet  temperature  was  (Table  1)  1000°  C with  a 
billet  velocity  of  500  mm/sec.  As  can  be  seen  the  maximum  temperature  rise  is  about 
30°  C in  the  worked  zone  as  shown:  Figure  26  shows  a simulation  identical  to  that  of 
Figure  25,  except  that  the  speed  has  been  doubled,  to  1000  mm/sec.  Again  the  maxi- 
mum temperature  rise  is  about  30°  C.  However,  as  can  be  noted  by  the  elongated 
1030°  C isotherm,  the  higher  temperature  zone  has  extended  back  into  the  die  and  for- 
ward into  the  product. 

(For  convenience.  Figure  27  taken  from  the  Battelle  report  is  presented.  The 
processing  conditions  are  identical  to  that  of  Figure  26,  a punch  velocity  of  1000  mm/s. 
The  only  difference  is  that  the  Battelle  run  was  terminated  at  . 145  sec  after  drawing 
was  initiated,  and  only  about  50%  of  the  product  had  been  formed.  Figure  was  run 
until  about  . 3 sec,  when  all  material  had  been  processed.  The  temperature  distribu- 
tions are  nearly  identical.  It  was  concluded  that  running  only  half  the  material  through 
the  die,  and,  hence  half  the  computer  expense,  was  sufficient  for  practical  purposes). 

Figure  28  shows  the  temperature  increase  starting  with  a cold  draw  (20°  C).  The 
punch  head  velocity  was  200  mm/s.  The  maximum  temperature  rise  is  about  320°  C, 
near  the  punch-shell  interface. 

Figure  29  shows  the  temperature  distribution  in  a cold,  40%,  draw  operation  (20° C 
ambient)  at  200  mm/s  using  an  optimal  double  curvature  die.  As  can  be  noted  by  com- 
paring Figures  25  and  30,  the  maximum  temperature  rise  is  higher  (by  about  50°  C)  in 
the  curved  die. 
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Figure  25.  Temperature  Distribution  in  Hot  Drawing  of  Steel  Shell  Through  a Conical  Die. 

(Billet  Temperature  1000°  C,  Tooling  Temperature  20°  C).  Velocity:  500  mm/s. 


Figure  26.  Temperature  Distribution  in  Hot  Drawing  of  Steel  Shell  Through  a Conical  Die. 

(Billet  Temperature  1000° C,  Tooling  Temperature  20°  C).  Velocity:  1000  mm/s. 


Temperature  Distributions  in  Hot  Drawing  of  Steel  Shells  Through  a Conical  Die 
0. 145  Seconds  After  Drawing  Starts 

(Billet  Temperature  = 1000°C,  Tooling  Temperature  ^ 20°  C)  Speed:  mm/sec 


28.  Temperature  Distribution  in  Drawing  of  Steel  Billet  Through  Various  Dies. 
(Billet.  Tooling  Initial  Temperature  20° C).  Velocity:  200  mm/s. 


Figure  29.  Temperature  Distribution  in  Drawing  of  Steel  Billet  Through  Various  Dies. 
(Billet.  Tooling  Initial  Temperature  20° C).  Speed:  200  mm/sec. 


Figure  3<T  Temperature  Distributions  in  Hot  Drawing  of  Steel  Shells  Through  a Double 
Curvature  Die  0.  145  Seconds  After  Drawing  Starts 

TUllet  Temperature  1000' C.  Tooling  Temperature  20°  C).  Speed:  mm/sec 


Figure  20  shows  the  hot  drawing  of  the  steel  billet  through  optimal  double  curvature 
dies.  The  conditions  used  were  identical  to  those  of  Figure  26  using  a conventional  die. 
As  can  be  noted  the  temperature  distributions  are  nearly  the  same,  they  differ  only  by 
about  10°  C. 

The  data  shown  in  Table  7 have  been  taken  from  the  paper  by  T.  Altan  and  F.  W. 
Boulger  in  Trans.  ASME,  J.  Engineering  for  Industry,  Vol.  95,  November  1975, 
pp  1009-1019.  It  should  be  noted  that  the  values  of  C and  m in  the  formula  for  hot 
working 


<T  = C (f)m 


are  functions  of  the  strain  and  temperature.  The  values  of  K and  n in  the  formula  for 
cold  working 

o = K (f)n 

are  functions  of  the  strain  rate. 

The  data  are  stored  in  memory  in  tabular  form.  The  computer  automatically 
choses  the  correct  values  of  C,  m,  K and  n by  interpolating  within  the  table. 
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In  this  study  mathematical  models  were  developed  and  used  to  establish  theoretical 
data  for  the  design  and  optimization  of  single  pass  die  rings.  Two  die  rings  were  con- 
sidered: (1)  optimal  curved  dies, and  (2)  conventional  tapered  rings.  The  effects  of 
temperature,  punch  speed,  friction,  heat  transfer,  and  material  were  used  in  develop- 
ing the  data  presented  in  the  section  entitled  Simulation  Results.  The  major  conclu- 
sions and  design  guidelines  for  each  of  the  two  die  types  are  given  below.  For  clarity 
also  given,  were  needed,  is  a specific  figure  or  table  number  from  the  Simulation  Re- 
sults section.  Also  undertaken  was  a study  to  show  the  effect  of  heat  transport  on  the 
drawing  process.  Conclusions  are  presented  at  the  end  of  this  section. 

I.  OPTIMAL  CURVED  DIES 


1.  The  design  of  optimal,  curved  die  profiles  has  been  demonstrated  theoreti- 
cally. As  shown  in  Figure  1,  these  dies  are  sigmoidal,  that  is  doubly  curved 
and  streamlined.  These  profiles  are  formed  from  the  blended  radii  of  two 
circles.  Optimal  in  the  present  context  means  that:  for  a specified  reduction, 
shell  metal  and  friction  condition,  this  curved  die  profile,  when  compared  to 

a conventional  tapered  die,  is  superior  in  that: 

a.  Less  energy  is  required  from  the  press  to  form  the  drawn  product. 

b.  The  residual  strain  distribution  in  the  product  wall,  at  the  die  exit,  is 
more  uniform  through  the  wall  thickness. 

The  first  factor,  less  energy  used,  indicates  a potential  cost  savings  b>  re- 
ducing plant  power  requirements.  While  the  findings  here  are  theoretical, 
an  experimental  study  conducted  by  US  Steel  for  drawing  sheet  and  strip,  in- 
dicated that  using  properly  designed  curved  dies  resulted  in  a 10%  energy 
savings.  The  second  factor,  homogeneous  strain  distribution  in  the  wall,  in- 
dicates a reduced  likelihood  of  failure  by  punch  through.  This  would  result 
in  reduced  scrap. 

2.  An  additional  consideration  in  the  die  design  is  the  potential  for  product  fail- 
ure due  to  punch  through  of  the  moving  mandrel.  In  this  study  a criterion  to 
predict  failure  was  postulated.  It  was  assumed  that  when  the  specific  draw' 
stress  exceeds  unity,  failure  occurred.  The  specific  draw  stress  was  found 
to  be  dependent  upon  that  temperaute,  friction,  mandrel  speed  and  the  stress- 
strain  curve  of  the  metal.  (See  Section  IV  on  Failure  Models  for  details  and 
limitations). 

3.  Computer  programs  which  perform  the  required  design  calculations  have  been 
developed.  The  programs  have  been  made  as  user-oriented  as  possible.  In- 
put is  punched  on  three  cards,  and  consists  oi  the  steel  being  used,  billet  size 
and  reduction  desired.  The  optimal  curved  die  profile  is  generated  in  about 
three  minutes  on  the  CDC  6500  at  a cost  of  about  $50.  The  program  auto- 
matically prints  out  the:  optimal  die  contour,  press  energy  requirements, 
loads  on  the  dies,  loads  on  the  shell,  loads  on  the  punch.  A failure  analysis 
and  diagnosis  of  potential  punch  through  is  automatically  performed. 
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4.  As  a result  ol'  this  study  it  was  possible  to  delineate  those  factors  which  do 

or  do  not  effect  the  optimal  die  shape.  It  was  found  that: 

a.  The  temperature  of  the  billet  has  no  effect  on  the  shape. 

b.  The  stress-strain  curve  of  the  billet  has  a small  effect  on  shape 
(Figure  2). 

c.  Friction  conditions  between  die  and  product  and  product  and  punch  have 
a small  effect  on  shape  (Figure  5). 

d.  The  velocity  of  the  punch  does  not  effect  shape. 

e.  The  most  significant  factor  effecting  shape  is  the  areal  reduction 
(Figures  6-14). 

5.  A study  to  develop  optimal  die  contours  in  drawing  a 105  mm  1045  steel  billet 

was  performed.  Figures  6-14  shown  are  the  contours  for  reductions  of  29%, 

32%,  36%,  40%,  41%Jt  42%,,  44%),  48%;  and  51%,  respectively. 

6.  It  was  possible  to  delineate  those  factors  which  effect  failure  in  optimal  dies. 

It  was  found  that: 

a.  Temperature  has  a significant  effect  on  raising  the  specific  draw  stress, 
which  determines  failure  (Figure  2,  Table  2).  It  was  found  that,  at  larger 
reductions,  certain  steels  cannot  sustain  even  a moderately  hot  draw. 

1095  and  52100  steels  at  750°  C will  fail  in  a 40%;  draw.  These  same 
steels  will,  in  the  cold  (20DC)  condition,  probably  sustain  a cold  draw 
without  failure. 

b.  Friction  has  a large  effect  on  potential  failure  (Figure  3 and  4).  It  was 
established  that  reducing  the  friction  between  shell  and  die  permits  larger 
reductions.  It  was  found,  conversely,  that  increasing  the  friction  between 
die  and  mandrel  allows  a larger  reduction.  This  latter  effect  may  be  ex- 
plained by  noting  that  the  force  necessary  to  draw  a shell  is  the  sum  of 
two  parts:  (a)  the  draw  stress  in  the  shell  wall,  caused  by  the  action  of 
the  mandrel  on  the  bottom  of  the  shell  and  (b)  the  friction  between  the 
mandrel  and  shell.  Thus,  when  the  friction  is  increased,  a greater  part 
of  the  total  load  is  carried  by  the  mandrel,  and  the  draw  stress  is  lowered. 

It  was  also  established  that,  regardless  of  friction,  40-41%  is  the  largest 
possible  reduction  in  a single  pass,  cold  draw  of  a 105  mm  projectile. 
Reductions  beyond  this  point  led  to  a specific  draw  stress  much  greater 
than  one.  It  was  also  shown  that  the  greater  the  reduction,  the  lower 
the  strain  homogeneity  in  the  product  wall. 
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n.  CONVENTIONAL  TAPERED  DIES 

1.  Unlike  the  optimal  curved  dies,  in  which  the  shape  is  not  known  a priority  but 
is  determined  as  part  of  the  problem,  conventional  dies  have  a fixed  known 
geometry  (Figure  A-l).  Thus,  the  computer  program  developed  to  accommo- 
date the  calculations  was  much  simpler.  Using  this  program,  the  effects  of 
friction,  temperature  and  mandrel  speed  were  delineated. 

2.  In  cold  drawing  regardless  of  the  friction  conditions  between  shell  and  die,  an 
increasing  die  angle  reduces  the  specific  draw  stress  and  hence  avoids  failure. 
Decreasing  the  shell-die  friction  enhances  the  process  by  lowering  the  likli- 
hood  of  failure  (Figure  15). 

3.  For  any  friction  conditions  in  cold  drawing,  there  is  a minimum  angle  below 
which  drawing  is  not  possible  (Figures  15  and  16). 

4.  Increasing  the  mandrel-shell  friction  in  cold  drawing  decreases  the  specific 
draw  stress  and  thus  lowers  the  chance  of  failure.  It  should  be  noted  that  the 
same  conclusion  was  drawn  for  curved  dies,  and  has  been  discussed  above. 

5.  Provided  the  two  coefficients  of  friction  are  equal,  hot  drawing  (up  to  1000°  C) 
is  very  similar  to  cold  drawing.  In  this  condition  the  punch  velocity  has  no 
effect  on  drawability.  The  condition  of  u-^  = is  not  met,  however,  in 
practice. 

6.  When  the  coefficients  of  friction  differ,  even  by  a small  amount,  the  draw 
stress  is  considerably  reduced  in  hot  drawing.  An  increased  punch  velocity 
will,  it  is  anticipated,  increase  the  liklihood  of  punch  through. 

7.  As  in  cold  drawing,  increasing  the  die-shell  friction  lowers  the  draw  stress. 

III.  HEAT  TRANSPORT 


A study  to  delineate  the  effects  of  heat  transfer  during  the  drawing  of  both  optimal 
and  conventional  dies  was  performed.  Both  hot  and  initially  cold  billets  were  investi- 
gated. The  following  conclusions  were  reached: 

1.  In  hot  drawing  (billet  uniform  initial  temperature  1000°  C)  the  temperature  in- 
crease, due  to  the  conversion  of  plastic  work  to  thermal  energy,  is  small,  on 
the  order  of  30°  C.  The  optimal  curved  die  has  a slightly  lower  maximum  final 
temperature  (1020°  C)  than  the  conventional  die  1030°).  It  was  concluded  that 
this  increase  is  insignificant  in  determining  the  product's  degraded  yield 
strength,  and  is  not  important  to  failure.  (Figures  25-27  and  30). 

2.  In  cold  drawing  (billet  uniform  initial  temperature  20°  C)  the  temperature  in- 
creases is  higher  than  in  hot  drawing  - it  is  between  330°  and  380°  C.  The 
curved  die  achieves  a higher  final  temperature,  400°  C,  than  does  the  tapered 
die,  350°  C.  These  temperatures,  at  normal  operating  speeds,  are  not  signi- 
ficant to  failure. 
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With  these  facts  in  mind,  the  broad  contributions  which  can  result  from  math 
modeling  studies  include: 

1.  A method  to  predict  the  suitability  of  present  equipment  to  future  shell  metals 
or  configurations. 

2.  A method  to  predict  the  effects  on  present  equipment  by  future  changes  in 
speed,  temperature  and  lubricants, 

3.  An  automated  method  . <•  the  preparation  of  technical  data  packages,  and  an 
aid  in  evaluating  vendor's  proposed  draw  presses  and  tooling  configurations. 

4.  A reduction  in  scrap,  and  thus  a subsequent  increased  economy  of  operation. 

Those  specific  processing  problem  areas  or  data  gaps  which  can  be  addressed 
by  the  results  of  this  task  include: 

1.  The  effects  of  drawing  of  the  higher  mandrel  speeds  from  a mechanical  press. 

2.  The  effects  of  temperature  in  a hot  draw  or  reduced  mechanical  properties 
in  the  product  wall. 

3.  The  effects  of  friction  between  the  punch  and  wall  and  between  the  wall  and 
die  rings. 

RECOMMENDATIONS 

Based  on  the  results  of  this  work,  it  is  recommended  that: 

1.  An  experimental  study  by  undertaken  to  confirm  the  optimality  of  the  sigmoidal 
die  rings  theoretically  designed  in  this  report.  Subject  to  such  confirmation, 
these  die  rings  should  be  made  available  for  contractor  use  in  shell  manufacture. 

2.  The  design  procedures  and  associated  software  for  shell  drawing  be  made 
available  to  engineers  involved  in  manufacturing  process  design. 

3.  The  design  procedures  be  extended  to  include  tandem  drawing  with  multiple 
die  rings  and  a tapered  mandrel. 

4.  The  failure  model  will  be  modified  to  include  plastic  instability. 

5.  Additional  work  on  hot  drawing  with  conventional  tapered  dies  with  unequal 
friction  coefficients  be  undertaken. 

(j.  New  modeling  studies  covering  each  step  in  shell  manufacture  be  undertaken. 
These  studies  should  include:  nosing,  billet  separation,  sizing,  cabbage  and 
pierce. 
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APPENDIX  A 

MATHEMATICAL  MODELING 
NOMENCLATURE 

a axial  stress 

p pressure 

A "(r|  - r^),  crossectional  shell  area 

r^  die  radius,  slab  method 

T2  die  radius,  slab  method 

x axial  coordinate  for  generalized  draw 

coefficient  friction  between  shell  and  mandrel 
^2*  coefficient  of  friction  between  shell  and  die 
P(x)  term 

Q(x)  term 

C constant 

Rq  radius 

tQ  thickness 

t^  thickness 

rQ  radius 

a die  semi-angle 

B constant 

Rj  radius,  optimal  curved  die 

R2  radius,  optimal  curved  due 

f effective  plastic  strain 

\ effective  plastic  strain  rate 

Z axial  coordinate  for  a tapered  die 

h thickness  for  a tapered  die 

Za  axial  position 

♦Explanation  next  page 
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Z.  axial  position 

t time 

nij  friction  factor  between  die  and  shell 

m2  friction  factor  between  shell  and  mandrel 

*In  the  computer  program  SHELL  friction  is  accounted  for  through  the  assumption  of 
Coulomb's  Law, 


tangential  shear  stress  _ T 
normal  pressure  p 


v 


where  is  the  coefficient  of  friction,  which  is  taken  as  constant  for  a given  die  and 
shell. 


In  the  programs  CDVEL,  SDVEL  and  SDTEM  another  expression  for  friction  is 
used.  It  is  assumed  that 

tangential  shear  stress 

— — yield  strength  of  the  shell  J y~3~ 

VT 

where  m is  the  friction  factor.  This  assumption  is  related  to  the  von  Mises  yield 
criterion. 


I.  BACKGROUND/RATIONALE 


Traditionally,  the  hot  or  cold  drawing  of  shell  has  been  carried  out  in  hydraulic 
presses  at  moderate  speeds  up  to  300  inches/min.  Recently,  presses  which  can  draw 
at  speeds  beyond  2000  inches/min  have  been  under  consideration.  The  die  and  process 
design  experience  is  largely  regulated  to  lower  speed  machines.  When  new  processing 
techniques  and  higher  speeds  are  introduced,  previous  tooling  experience  alone  is  not 
sufficient  to  design  the  drawing  process  adequately.  For  example,  when  the  drawing 
speed  is  increased  5 to  10  times  over  conventional  speeds,  it  is  necessary  to  utilize 
quantitative  engineering  analysis  in  order  to  predict  under  which  conditions  failure 
such  as  punch-through  or  wall-tearing  could  occur.  This  problem  is  typical  of  the 
interaction  of  temperature  and  strain  rate  in  reducing  strength  at  high  processing 
speeds.  Two  objectives  of  this  task,  then  are  to  provide  models  which  can  predict 
potential  failure,  and  establish  optimum  combinations  of  the  processing  variables. 

The  mathematical  techniques  required  to  meet  these  objectives  are  discussed  below. 

A review  of  the  technical  literature  on  the  modeling  of  the  shell  drawing  process 
indicates  that,  basically,  the  following  three  methods  of  attack  have  been  applied  to 
the  mathematical  study  of  the  process:  (a)  the  "Slab"  or  "Sachs"  Method,  (b)  the 
"Method  of  Characteristics"  or  the  "Slipline  Analysis",  and  (c)  "The  Upper  Bound"  or 
"Minimum  Energy"  Method. 

Using  the  slab  method,  Snchs(2)  and  his  coworkers  conducted  the  first  mathe- 
matical analysis  of  the  shell  drawing  process.  Their  study  treated  the  process  as  a 
plane-strain  problem  and  it  was  based  on  three  main  assumptions:  (1)  the  deformation 
was  a simple  uniform  (uni-axial)  extension,  (2)  the  principal  stress  directions  were 
everywhere  perpendicular  and  parallel  to  the  punch  surface,  and  (3)  the  plastic  region 
extended  between  the  two  transverse  planes  containing  the  circles  of  contact  of  shell 
and  die  at  entry  and  exit.  W'ith  these  assumptions,  Sachs  et  al.  , (2)  were  able  to  de- 
rive and  graphically  solve  the  equilibrium  equations  for  straight  as  well  as  curved  dies. 
He  also  crudely  introduced  the  effect  of  strain-hardening  and  of  the  strain-rate  depen- 
dency of  the  shell  material.  Thus,  Sachs'  plane-strain  analysis  allows  one  to  deter- 
mine, within  these  approximations,  the  effect  of  process  variables  upon  the  draw  force 
and  upon  the  drawing  limits  in  drawing  very  thin-walled  shells.  However , this  analysis 
is  not  adequate  for  drawing  of  medium  or  thickwalled  shells  (artillery  shells)  and  it 
does  not  include  the  thermal  effects,  i.  e. , heat  generation  and  heat  transfer,  which 
are  important  aspects  of  the  deformation  process  at  higher  speeds. 

The  slipline  analysis  assumes  a rigid-perfectly  plastic  material.  This  method  has 
been  originally  applied  by  Hill^  to  the  solution  of  plane-strain  shell  drawing  through 
conical  dies.  Similar  studies  have  also  been  conducted  for  curved  dies  by  Green, 
Richmond  and  Devenpeck^>  and  Sowerby,  et  al.  , (®).  Recently,  the  slipline  analy- 
sis has  been  extended  to  study  axisymmetrio  extrusion  and  drawing  operations  by 
Richmond.  Although  the  slipline  analysis  is  among  the  most  rigorous  methods  for  in- 
vestigating metal  forming  problems,  it  essentially  neglects  strain-hardening  effects, 
strain-rate  effects,  as  well  as  thermal  phenomena  in  the  deforming  material.  Con- 
sequently, its  application  to  optimizing  actual  shell  and  cartridge  drawing  operations 


requires  modification.  The  outstanding  advantage  of  this  method  is  that  the  stress  field 
is  very  accurately  determined.  Failures,  such  as  cracking  or  central  burst,  depend 
on  an  accurate  estimate  of  the  state  of  stress  for  their  prediction.  Thus,  this  tech- 
nique is  required  for  examining  failure  modes.  Progress  utilizing  this  technique  is  pre- 
sented in  this  report. 

The  upper-bound  method  describes  the  metal  flow  in  the  shell  drawing  operation  by 
considering  an  admissible  velocity  field  that  satisfies  the  incompressibility,  continuity, 
and  the  velocity  boundary  conditions.  Based  on  this  velocity  field,  the  deformation, 
the  shear  (if  velocity  discontinuities  are  present),  and  the  friction  energies  are  com- 
puted to  give  the  total  forming  energy  and  also  the  forming  load.  Based  on  limit  the- 
orems, (10)  this  calculated  forming  load  represents  an  upper-bound  to  the  actual  form- 
ing load.  Often,  the  velocity  field  considered  includes  one  or  more  free  parameters 
that  are  to  be  determined  by  minimizing  the  total  forming  energy  with  respect  to  those 
parameters.  Thus,  the  determined  values  of  the  parameters  enable  abetter  upper- 
bound  velocity  field  to  be  found. 

Most  of  the  upper-bound  solutions  available  for  processing  of  tubular  products  have 
been  developed  for  extrusion.  However,  because  of  geometrical  and  fundamental  sim- 
ilarities in  metal  flow,  these  analyses  can  also  be  applied,  with  appropriate  modifica- 
tions, to  the  shell-drawing  process.  A radial  velocity  field  for  tube  extrusion  through 
conical  dies  was  derived  by  Avitzur(H)  anf]  aiso  applied  1o  practical  extrusion  studies 
by  DePierre.  (1^)  Chen  and  Ling^)  also  utilized  the  upper-bound  approach  extensively 
to  obtain  solutions  to  axisymmetric  extrusion  problems  with  various  curved  dies.  Chang 
and  Choi^4)  conducted  similar  studies  on  the  extrusion  of  tubes  through  conical  and 
curved  flies.  Mehta,  et  al.,  0®)  developed  a unique  superposition  technique  to  develop 
a velocity  field  without  discontinuities  in  analyzing  the  tube  extrusion  through  conical 
dies. 


The  upper-bound  solutions,  discussed  above,  assumed  that  the  material  was  rigid- 
perfectly  plastic,  ft  is  possible  to  introduce,  with  some  difficulty,  the  strain-hardening 
and  strain-rate  effects  into  these  solutions.  However,  with  the  exception  of  Mehta's 
solution,  all  other  upper-bound  analyses  include  discontinuities  in  their  velocity  fields. 
Therefore,  these  discontinuities  result  in  artificial  and  unrealistic  sources  of  heat 
generation  when  a thermal  analysis  of  the  drawing  process  is  conducted.  Consequent- 
ly, the  existing  upper-bound  studies  required  significant  improvement  in  order  to  de- 
velop a sound  mathematical  model  of  the  shell  drawing  process.  That  modification  is 
given  in  this  report  and  in  greater  detail  in  the  Battelle  report. 

Although  the  influence  of  temperature  and  speed  in  shell  drawing  processes  is  well 
known,  a realistic  thermal  analysis  of  the  shell  drawing  operation  does  not  exist. 
However,  several  attempts  have  been  made  in  the  past  to  predict,  at  least  approx- 
imately, the  temperature  distributions  in  some  kinematically  steady-state  deforma- 
tion processes,  such  as  axisymmetric  extrusion  and  drawing  of  rod  and  wire.  Most 
recently,  with  the  wide  use  of  digital  computers,  numerical  methods  have  been  more 
widely  introduced.  These  methods  are  capable  of  describing  the  complex  physical 
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phenomena  involved  in  extrusion  or  drawing  without  the  many  simplifying  assump- 
tions as  it  is  the  case  in  using  analytical  methods  of  solution. 

in  plane-strain  extrusion,  Bishop(16)  developed  the  principles  of  a numerical 
method  by  separating  the  complex  problem  of  heat  generation,  conduction  and  trans- 
portation into  two  parts.  Thus,  the  heat  generation  and  transportation  are  regarded 
as  occurring  instantaneously,  followed  by  an  interval  in  which  conduction  takes  place 
as  for  a stationary  medium.  Altan  and  Kobayashi^17)  developed  this  method  further 
and  developed  computer  programs  for  determining  temperature  distributions  in  ex- 
trusion of  round  bars  by  using  experimentally  obtained  data  on  metal  flow. 

Altan' s(17>  18>  19»  20)  studies  have  been  further  expanded  by  using  an  axisymmetric 
upper-bound  velocity  field  without  discontinuities  to  calculate  the  heat  generation  and 
the  heat  transport  during  deformation.  This  work  has  been  further  applied  to  pre- 
dicting temperatures  in  extrusion  of  rods  and  in  drawing  of  rod  and  wire. 

Based  on  the  relative  advantages  and  disadvantages  of  the  available  methods  dis- 
cussed above,  it  was  decided  to  simultaneously: 

1.  Develop  a very  general  upper-bound  approach  which  could  handle  high  and  low 
forming  rates  and  temperature  generation. 

2.  Expand  the  slip  line  (method  of  characteristics)  techniques  to  include  high 
forming  rates. 

3.  Expand  the  Sachs'  or  Slab  Method  for  failure  prediction. 

The  first  effort  was  carried  out  under  contract  with  Battelle.  The  second  was  per 
formed  in-house  at  Frankford  Arsenal.  The  third  was  a joint  Battelle/Frankford  Arse 
nal  effort. 


II.  THEORY 


In  order  to  make  the  present  report  complete  and  self-contained,  the  basic  theo- 
retical considerations  of  all  methods  used  in  developing  the  various  math  models  are 
sketched  here.  For  details  concerning  the  upper  bound  method,  together  with  the  as- 
sociated computer  programs  reference  to  the  Battelle  report  is  made. 

IT  A.  THE  UPPER  BOUND  METHOD 

T he  basis  of  the  method,  due  originally  to  Prager^7),  and  modified  by  Altan (*7) 
and  his  associates  is  essentially  as  follows: 

1 A velocity  for  the  flowing  metal  is  guessed.  The  "best"  guess  is  given  by  ad- 
hering to  the  following  principals: 

a.  The  flow  must  be  incompressible. 

b.  The  flow  must  be  continuous. 

c.  The  flow  must  satisfy  the  velocity  boundary  conditions. 
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2.  Based  on  the  velocity  field  generated  above,  the  deformation  (strain  rates, 
velocity  gradients,  etc  ),  the  shear  energy  (if  velocity  discontinuities  were 
allowed),  and  friction  energies  are  calculated.  Based  on  the  limit  theorems 
of  plasticity^17)  the  forming  load  is  higher  than  the  actual  load  and  it,  there- 
by, represents  an  upper  bound  to  the  actual  forming  load.  It  also  follows 
that  the  lower  the  upper  bound  load,  the  better  prediction  of  metal  flow.  One 
way  of  insuring  a "best"  (smallest)  upper  bound  is  to  use  the  calculus  of  var- 
iations. Unknown  parameters  are  introduced  into  the  energy  function  and 
the  entire  function  minimized  according  to  the  rules  of  the  calculus  of  var- 
iations . 

As  an  example  and  to  develop  the  theory,  the  case  of  drawing  through  conical  dies 
is  presented.  As  shown  in  Figure  A-l  a tubular  billet  of  outer  radius  1 + r0  and  an 
inner  radius,  rQ,  is  drawn  through  a conical  die  of  half  angle  a,  to  a final  wall  thick- 
ness b,  using  a cylindrical  punch  of  radius  rQ. 


A basic  velocity  field  is  assumed  by  a plastically  deforming  region  (region  2 in 
Figure  A-l)  bounded  by  the  velocity  discontinuity  surfaces  ACQ  and  BDP.  The  flow 
lines  in  the  regions  before  and  after  the  deformation  region  are  also  assumed  to  be 
straight  lines,  parallel  to  the  axis  of  drawing.  The  radial  distances  R|<  and  r^,  of  the 
flow  lines  from  the  surface  of  the  punch  are  related  to  each  other  because  of  incom- 
pressibility. Thus: 
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where  Ve  is  the  velocity  at  the  entrance  or  billet  velocity.  Unlike  tube  extrusion,  in 
shell  drawing,  the  process  speed  is  controlled  by  the  punch  velocity,  which  is  the 
same  as  the  product  velocity.  The  billet  velocity,  Ve  for  unit  punch  velocity,  Vf,  is 
given  by 
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be  represented  by  the  equation 
G rk  -•  II  , (3) 

where  E,  F,  G,  II  and  0 are  to  be  determined  from  the  boundary  conditions.  Omitting 
the  details,  which  are  given  in  Reference  1.  The  net  result  is 


Let  the  flow  lines  in  region  2,  Figure  A-l, 
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Figure  A-l.  Geometry  of  Axisymmetric  Metal  Deformation  in  Shell  Drawing 


■ 


in  region  (1),  r = - r 
k 


in  region  (2),  r = 


o'V'oH2-' o2) 


V , 
e 


1 


k z - r cot~a 

r - r b ’ 


o 


and  in  region  (3),  r = r - r 

K O 


(4  . a) 
(4  b) 

(4  c) 


The  intersection  between  the  flow  lines  in  region  (2)  and  region  (1)  determines  the 
equation  for  the  line  ACQ  of  velocity  discontinuity.  This  equation  is  given  by 
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The  parameter  dj  in  Equation  4,  separating  the  flow  regions,  is  not  as  yet  determined 
by  imposition  at  the  boundary  conditions.  In  the  spirit  of  the  calculus  of  variations  the 
parameter  dj  will  be  determined  by  the  condition  that  a certain  integral  be  a minimum. 

For  a given  extrusion  ratio,  die  angle,  and  friction  factor,  the  actual  external 
work  rate  is  a function  of  the  only  unknown  parameter  dj  of  Figure  A-l.  The  unknown 
parameter  d^  can  be  derived  for  the  determination  of  an  optimum  velocity  field,  using 
the  upper-bound  theory.  For  this  purpose,  dj  must  be  such  that  the  rate  of  energy  E, 

The  energy  rate,  E,  is  given  by 

/f  • AVdS  - ft  -AVdsl,  (5) 
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consumed  in  the  process,  is  a minimum. 
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flow  stress  of  shell  material 
volume  of  the  deformation  zone 

friction  factor  at  the  shell-die  interface,  (0  < f ^ < 0.577) 
velocity  difference  at  the  shell-die  interface 
shell-die  interface  area 

friction  factor  at  the  shell-punch  interface,  (0  ^ f ^ 0.577) 
velocity  difference  at  the  shell  punch  interface 
shell-punch  interface  aren  excluding  the  product 
strain  rates. 


fn  order  to  evaluate  Equation  5 we  need  to  substitute  expressions  for  the  velocities  and 
strain  rates. 
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For  incompressible  materials,  the  velocity  fields  in  steady-state  axisymmetric 
plastic  deformation  processes  can  be  defined  by  introducing  a flow  function  such 
that 

l '<)  \j/  , 1 <i  4/ 

U = — — - and  V — (6) 

r dz  r dr  w 

where  u and  v are  the  velocity  components  in  the  radial  (r  -)  and  the  axial  (z-)  direc- 
tions, respectively,  and  = 4>  (r,  z) . The  flow  lines  are  then  given  by, 

* = Ck  (7) 

where  is  a constant  along  the  k-th  flow  line  and  proportional  to  the  rate  of  volume 
flowing  through  the  flow  layer.  The  volume  rate  flowing  through  a flow  layer  bounded 
by  the  k-th  flow  line  is  proportional  to-£-Ve  r^  (r^  + 2 r0)  where  r^,  in  various  regions, 
is  given  by  Equation  1.  The  flow  function  of  the  basic  flow  pattern,  is  expressed 
in  the  form  of 

* = T (8) 

The  velocity  components  are  then  calculated  according  to 
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The  velocity  field  given  by  Equation  8 contains  velocity  discontinuities.  A continuous 
velocity  field  can  be  obtained  by  superimposing  several  basic  flow  patterns,  each  of 
which  is  given  by  a pair  of  velocity  discontinuity  lines  as  shown  in  Figure  A-l.  The 
details  are  given  in  Reference  1 . The  strain  rates  may  be  in  turn  calculated  from 
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Once  the  optimum  value  of  dj  corresponding  to  the  minimum  energy  rate,  E,  is 
known,  then  the  best  member  of  the  selected  class  of  velocity  fields  is  given  by  Equa- 
tion 9.  These  calculations,  including  the  minimization  of  E (dj),  have  been  program- 
med for  use  at  Frankford  Arsenal  or.  a Cl)C  9500/0000  computer 

II  B.  PERTURBATION  APPROACH 

The  objectives  of  this  portion  of  the  work  were  to  develop  mathematical  methods 
capable  of  providing:  (a)  a detailed  stress  analysis,  (b)  a failure  analysis  of  drawing, 
and  (c)  solutions  to  problems  not  handled  by  the  upper  bound  method 

In  order  to  meet  the  above  objectives  the  work  has  been  divided  into  three  sub- 
tasks. A brief  description  of  these  sub-tasks  follows. 

Sub-Task  1.  Evaluation  of  Existing  Perturbation  Methods 

This  sub-task  consisted  of  exploring  perturbation  techniques  which  have  proven 
successful  in  other  areas  of  deformation  mechanics.  The  possibility  of  application  of 
this  method  to  plasticity,  and  to  metalforming  problems  in  general,  was  found  to  be 
excellent.  Lin^^  and  Fox^^  developed  the  technique  for  compressible  fluid  mechanics. 
Richmond(23j  later  noted  the  similarity  between  the  fluid  mechanics  case  and  plasticity. 

It  was  apparent  that  this  technique  is  capable  of  giving  a detailed  stress  analysis  of  the 
deforming  product.  This  capability  is  important  in  determining  failure  modes  in  draw- 
ing. This,  of  course,  cannot  be  done  using  the  upper  bound  method. 

A second  perturbation  method,  suggested  by,  Onat(2i)  includes  the  effects  of  work- 
hardening and  flow  stability  on  plasticity.  A careful  examination  of  these  two  preoedures 
was  completed,  and  the  two  techniques  were  persued . Details  of  this  investigation  are 
presented  below. 

Sub-Task  2.  Adaptation  of  Lin's  Method  for  Metal-Forming  Process 

This  model  is  applicable  to  determining  residual  stress  and  strain  during  process- 
ing. The  model  is  idealized  in  that:  the  product  metal  is  cold,  non-work-hardening 
and  experiences  no  friction  between  the  die-shell  and  mandrel-shell  interfaces. 

Sub-Task  3.  Adaptation  of  Gnat's  Method  for  the  Shell  Forming  Process 

This  model  is  used  in  conjunction  with  Lin's  method,  to  study  the  states  of  stress, 
strain  and  velocity  in  the  product.  The  model  is  less  idealized  than  that  of  sub-Task  2 
in  that  work-hardening  can  be  considered. 

The  basis  of  the  Lin  method  is  presented  in  References  21,  22.  The  concept  is  as 
follows : 


1.  The  partial  differential  equations  governing  the  plastic  flow  are  written  down 
in  terms  of  the  characteristics  or  slip  lines.  If  the  equations  are  not  hyper- 
bolic the  method  does  not  apply. 


2.  The  independent  variables  are  the  characteristics  of  the  original  field  equa- 
tions. The  dependent  variables  are  the  stresses,  velocities  and  Eulerian 
positions. 

ese  dependent  variables  expanded  in  a series  involving  a small  (perturbation) 
parameter  e , and  each  term  in  the  series  is  treated  as  a function  of  the  character- 
istic parameters. 

In  order  to  illustrate  the  method,  the  problem  of  tubular  drawing  through  dies  is 
sketched  below.  The  notation  is  essentially  that  of  Reference  23. 

Under  the  assumption  that  the  deforming  material  is  rigid -perfectly  plastic,  the 


characteristic  equations  are 

(P  + 0 ),  q + (Z  ± r),a/2r  = 0 (11. a) 

(P  - 0 ),  p - (Z  + r),  p/2r  = 0 (11  .b) 

Z,  a = tan  0 r,  a (11.  c) 

Z,  (3  = cot  0 r,  (3  (11. d) 

U,  or  - V0,  a + (U  , a - V ,a)/2r  = 0 (11  .e) 

r z 

V,  0 + U0,  0 + (U  , 0 + V ,0)/2r  = 0 (11.  f) 

z r 


In  Equation  11:  r and  z are  cylindrical  coordinates,  a and  p are  the  characteristic 
parameters,  and  P,  0,  u and  v are  related  to  the  usual  stress  and  velocity  components 


by 

p = - (or  + „z)/4k  (12. a) 

20  = tan 1 [ (0  - 0 )/2r  ] (12. b) 

z r rz 

U = u cos  0 + v sin  0 (12.  c) 

V = - u sin  0 + v cos  0 (12. d) 

The  following  assumptions  are  now  made: 
o 2 

P = P + e P’  +6  P"  + . . . (13. a) 

0 = 0°  +e0'  + e^0"  + ...  (13.  b) 

o 2 

r = r + e r'  + e Z"  + . . . (13. c) 

z = z°  + e z'  + e ^ Z"  + . . . (13. d) 
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where  each  of  the  P,  0 , r,  7 are  functions  of  the  true  characteristics  a and  0. 


Equation  13  is  now  introduced  into  Equation  11  and  like  order  terms  collected.  The 
zeroth  order  perturbed  counterpart  of  Equation  11  a-d  is  thus: 


00  .00 

(cos  0 ) Z , a = (Sin  0 ) r , a 

(14.  a) 

(P°  + 0 °),  a = 0 

(14. b) 

0 0 

(p  - 0 ),  a = 0 

(14.  c) 

00  00. 

(sin  0 ) Z , p - - (cos  0 ) r , P 

(14.  d) 

The  first  and  higher  order  perturbed  quantities  have  not  as  yet  been  derived. 

By  examining  Equation  14,  it  can  be  seen  that  they  represent  a simpler  system  of 
equations  than  does  Equation  11.  While  Equation  14  is  not  linear  the  terms  involving 
1/r  from  Equation  11  do  not  appear.  It  is  anticipated  that  Equation  14  and  the  higher 
order  perturbed  equations  will  lead  to  a solvable  system  of  equations. 

II  C.  THE  ONAT  PERTURBATION  METHOD 

The  basis  of  the  method  is  presented  in  Reference  (24).  A derivation  of  the  funda- 
mental equations  is  presented  below. 

For  a rigid-work  hardening  body  the  stress-strain  relations  are 

;i)  " FS...i2,  ,i2>0  (15. a) 

‘ij  ° , if  J2<  k2,  J2<  0 (15. b) 

where 


stress  tensor 

strain  rate  tensor 

stress  deviator  = 0 - 1/3  „ y 

^ PP  ij 

1/2  S..  S.. 
ij  iJ 

material  derivative  of  J9 
yield  stress  in  shear 


It  is  assumed  that  the  velocities,  stress  and  strain  rates,  referred  to  convention- 
al Cartesian  coordinates,  Xj,  can  be  expanded  in  a time  series  as 


tr-f- 

II 

v;  <v 

+ t 7.  (a  ) 
i ' p7 

+. . . + 

(16.  a) 

••  (a  , t)  = 
oi]  p 

°ij  (ap) 

+ t o'.  (a  ) 
iJ  p' 

+. . .+ 

(16.  b) 

..  (a  , t) 
P 

e°ij  <V 

+ t e'..  (a  ) 
i]  p7 

+.  . . + 

(16.  c) 

where 

t i = Lagrangian  acceleration 


°ij  = rate  of  stress 

e'jj  = velocity  strain  (strain  rate) 

a^  = Lagrangian  coordinates 

and  the  quantities  with  a superscript  "0"  are  evaluated  at  t = 0,  the  yield  point  state. 
These  satisfy  the  following  equilibrium,  yield  condition  and  kinematic  quantities 


_ o 

O ..  , . = 0 

1J  J 


1/2  S.  ° S.  ° = k2 
i]  iJ 


1/2  ^.J0  + Vj,i°>  ^ F^Sij°  ^ 


The  first  order  quantities  satisfy 


f = 1/2  (7.  .+>..)-  1/2  (V.  ° V °+V.  ° V ) 

D bJ  J*1  bP  PfJ  ],P  p,i' 


o'..  . - V ° n ..  ° = 0 

b.J  P,J  iJ.P 


(17) 

(18) 
(19) 


(20) 


Equations  18,  19,  29  and  22  govern  the  first  order  quantities. 


A comparison  of  Upper  Bound  and  Perturbation  methods  is  shown  in  Tables  A-l 
and  A-2. 

f 

j Table  A-l . The  Range  of  Application  for  Each  Math  Model 

r 

Metalworking  Application 


Math  Model 

Stress 

Analysis 

Forming 

Loads 

Failure 

(Punch  Through) 

Optimum 
Die  Design 

Upper  Bound 

No 

Yes 

Limited 

Yes 

Lin  Perturbation 

Yes 

Limited 

Limited 

No 

Onat  Perturbation 

Yes 

Limited 

Yes 

No 

Table  A-2.  The  Range  of  Metalworking  Parameters  for  Each  Math  Model 

Metalworking  Parameter 


Math  Model 

Elevated 

Temperature 

Strain 

Hardening 

Strain 

Rate 

Friction* 

Upper  Bound 

Yes 

Yes 

Yes 

Yes 

Lin  Perturbation 

No 

No 

No 

Yes 

Onat  Perturbation 

No 

Yes 

No 

?** 

* Between  tooling  and  workpiece. 

‘+  Theory  is  insufficiently  developed  to  comment. 
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II.  D SLAB  METHOD 

In  the  preceding  sections,  the  upper  bound  method  has  been  used  to  analyze  draw- 
ing. As  mentioned  there,  the  upper  bound  method  suffers  from  the  disadvantage  that 
the  stress  distributions  calculated  are  not  accurate.  Failure,  if  it  occurs,  is  due  to 
tearing  of  the  shell  wall.  Tearing  in  turn  is  related  to  an  excessive  axial  stress  in  the 
drawn  product.  Thus,  for  purposes  of  design  it  is  necessary  to  have  an  accurate  es- 
timate of  the  stress  in  the  product  wall.  The  purpose  of  this  section  is  to  present  an 
analysis  which  meets  this  requirement. 

Since  it  has  been  established  earlier  in  this  report  that  optimal  dies  are  dies  which 
are  curved  rather  than  conical,  the  stress  analysis  formulae  given  here  are  very  gen- 
eral and  are  applicable  to  any  curved  or  conical  die.  Also,  since  in  practice  the 
mandrel  is  tapered  or  curved  this  also  is  considered. 

Figure  A-2  shows  a diagram  of  the  shell  drawing  process.  Both  the  die  contour 
and  mandrel  contour  are  allowed  to  be  specified  functions  of  the  axial  coordinate,  x. 
Since  the  shell  geometry  has  cylindrical  symmetry  about  the  x-axis,  only  the  two 
radii,  rj  and  r2,  are  required  to  completely  specify  the  shape.  The  forces  or  stres- 
ses acting  on  the  element  (exploded  view)  are:  the  stress,  0 A;  the  normal  pressure, 
p , and  the  frictional  forces  p^p  and  P2P-  Equating  these  forces  to  zero  results  in  the 
following  equation  of  equilibrium: 


°A  + _ d ( o A ) dx 

dx 


°A  - pr2  dr2  2"  t-  pr^  dr^  2" 


+ P2p  2 r2  dx  2n  ~ ^l  Pi  dx  2"  0 

Simplifying  and  noting  that 


(23) 


A " (r22  - rj2) 
it  follows  from  (23)  that 

d 

dx 


124 ) 


- 

■ 

"(r22 

n2> 

+ 2p 

r9  dr2 

dri 

- r.  1 

2 

dx 

1 

dx 

+ 2 p 


- r 


(25) 


1^1  + r2  ^2 


0 


1 ■ 

i 

: 


! 


I 
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According  to  Sachs  (2-a),  the  yield  condition  is 
2 


p t o 


(26) 


V3 


where  o is  the  yield  strength  of  the  metal.  Eliminating  p from  Equation  25  by  using 
26  gives 


dp 


P(x)°  = Q(x) 


(27) 


dx 


where 


P(x>— 


2 (*2  r2  - "l  rl 


(r22  - r*) 


Q(x)=-  2 


(r2  r2'  " rl  ri  > + ("2  r2  " "1  rl> 


(28) 


2 2 
<r2  " rl  ) 


Equation  27  is  a first  order  linear  differential  equation  with  variable  coefficients 
for  the  unknown  stress  o(x).  The  solution  of  this  equation  is 


o(x)  = C exp  ( 


dx) 


exp  ( - / pd? 

0 . 

exp  ( - / pdx)  • Qdx 

J 

_ 

L J 

(29) 


where  C is  a constant  of  integration.  This  constant  may  be  found  by  noting  that  at  X 0, 
at  the  entrance  to  the  die,  the  applied  force  is  zero.  Thus 


o (o)  = 0 

is  the  required  condition. 


(30) 


The  above,  Equation  24,  is  valid  only  for  a single  die  ring.  In  the  case  of  two  or 
more  rings  (tandem  drawing)  this  would  be  replaced  by  the  conditions  that:  (a)  the 
stress  is  zero  going  into  the  first  die  and  (b)  the  stress  must  be  continuous  upon  entry 
into  the  second  die.  This  work  has  not  as  yet  been  carried  out. 
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In  general,  although  the  solution  given  by  Equation  24  is  in  closed  form,  the  actual 
integration  is  difficult  analytically  except  under  certain  special  cases.  In  general, 
numerical  integration  may  be  required.  For  the  conical  die  the  integration  is  simple 
and  is  given  in  the  next  section. 


The  failure  criterion  is  given  by  Sachs  as  follows:  When  the  total  drawing  stress 
(total  drawing  load/shell  area  at  die  exit)  exceeds  2/^3  times  the  yield  strength  (as 
measured  at  the  exit),  then  tensile  wall  failure  occurs.  The  total  drawing  stress  is  the 
sum  of  all  the  axial  forces  acting  on  the  shell  and  on  the  mandrel  divided  by  shell  cross- 
sectional  area  at  the  exit.  This  is  not  to  be  confused  with  the  "actual  stress"  o , given 
in  the  above  formulas.  The  drawing  stress  is  actually  a fictious  quantity  which  is  use- 
ful in  defining  the  drawing  limits  or  failure  in  a proposed  process. 

’n  this  report  the  special  case  of  a conventional  conical  was  treated.  In  order  to 
accomplish  this  the  above  equations  are  reduced  to  the  equations  for  a conical  die  of 
semi-angle  « and  a straight  mandrel.  Thus  referring  to  Figure  A-3,  the  inner  radius, 
r,,  is  a constant,  rQ,  or 

ri  = rQ  (31) 

and  the  radius  ^ is  a linear  function  of  Z 

R1  “ R0  Ri  ' R0 

r2  = ( -)  Z + Rx  - ( ) Za  (32) 

Za  - zb  za  - zb 

Using  Equations  31  and  32  in  the  previous  equations  results  in  a stress 
distribution 

1 + B 2 

o(Z)  ( ) o 

B V3 

where 


1 - ( 


r2  ' r0 


B 


(33) 


R0  ' r0 


B = - u y)/tan  a 


At  the  exit  where  r2  Rj  the  actual  stress  is 


n 


1 + B 


( ) 1 

B 


“I  - 

( ) 

R0  * r0 


(34) 


(35) 


The  total  draw  force  consists  of  two  parts.  Fj  and  1'2:  where  F,  is  the  resultant 
of  the  tensile  stresses  acting  on  the  shell  at  the  exit,  and  1'2  is  the  force  on  the  mandrel 
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due  to  the  components  of  normal  pressure  and  friction  acting  in  the  axial  direction. 
From  Equation 


2a  1 + B [ / R - r \ _ 

F1  .(R^-ro2) ( ) 1-  ( -) 

yfj  B \n  -r  j 


F2  = 2,7r0/p2  J 


dz  = dh/tana 


P = -|r  o-o 


Ro  - ro 


2 2 


F2  = 2,rr0 


tana 


f \ 1 + B , h B 


Ri  - ro 


which  upon  simplification  is 


dr  (40) 


F2  2"r0 tx( 

VI  ^ 2 


^ 1 # w 2 

If  ^2  then  this  is  modified  to  read 


* " l0 
f2  2,1  ro  ti  In  ( ) 

VI  tana  li 

^1  u2=  ^ 


The  drawing  stress  is  obtained  by  summing  F^  and  lr2  and  dividing  by  the  exit 
area  ” (Rj2  - rQ2).  t^us  the  specific  draw  stress  is 


• - > V > - • ‘ 


^ 2 ^ 


Equations  43  and  1 1 arc  the  design  equations  tor  potential  failure.  The  failure 
prediction  procedure  is  as  follows: 

1 . Equation  13  or  44  is  evaluated  using  the  appropriate  geometri  : and  frictional 
factors.  The  value  of  o used  must  be  the  tensile  yield  stres:  at  the  exit 
under  whatever  temperature  and  strain  rates  are  experienced  by  the  shell  at 
the  exit. 

2.  For  a cold  draw  at  low  speeds  o can  be  taken  as  the  ordinary  sequence  of 
values  of  f from  a stress-strain  curve.  In  this  case  only  wo  -k  hardening  is 
important;  that  is 

a ~ o ( e ) (45) 

3.  For  a hot  draw  at  moderate  or  high  speeds,  the  effects  of  str;  in-rale  at  the 
temperature  of  the  draw  is  important;  that  is 


( 5) 
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ii  ( f . f ) 


4.  Both  cases  mentioned  above  have  been  automatically  covered  in  one  computer 
program  SHELL.  In  this  program  Equations  43  and  44  are  evaluated  using 
estimates  of  the  effects  of  strain  rate  and  strain  hardening. 


III.  HEAT  TRANSFER 

As  mentioned  previously,  the  effect  of  temperature  on  a drawing  process  is  pri- 
marily seen  in  the  influence  of  strain  rate  on  the  lowering  of  metal  strength  at  elevated 
speeds  and  temperatures.  As  a general  conclusion:  The  higher  the  mandrel  speed 
(or  the  higher  the  strain-rate),  the  more  critical  is  the  operating  temperature.  As 
the  temperature  is  raised  not  only  is  the  metal  strength  degraded  but  other  physical 
properties  change  as  well.  Among  these  are  the  density,  specific  heat  and  thermal 
conductivity.  As  is  well  known  when  these  quantities  depend  on  temperature,  the  net 
resulting  effect  on  the  governing  equations  in  nonlinearity.  In  addition,  since  the 
temperature  will  vary  from  point  to  point  in  the  workpiece  so  will  the  various  prop- 
erties and,  hence,  the  conduction  of  heat.  In  order  to  assess  this  later  effect,  con- 
ductive inhomogeneity,  the  problem  of  heat  conduction  in  a nonhomogeneous  medium 
was  studied. 

Consider  the  case  in  the  linear  theory  in  which  heat  transfer  diffusion  depends 
upon  a single  spacial  coordinate,  x,  and  time,  t.  Then  the  fundamental  conservation 
equation  to  be  solved  in  the  cartesian,  cylindrical  or  spherical  coordinate  systems 
is  (25,26) 


where  m has  the  value  0,  1 or  2 according  as  to  whether  cartesian,  cylindrical  or 
spherical  symmetry  is  being  considered.  If  heat  conduction  is  being  considered  K is 
the  conductivity,  p is  the  density  and  C the  specific  heat.  If  diffusion  is  being  con- 
sidered, the  produce  p C is  unity  and  K is  the  diffusion  constant.  K,  p and  C are 
assumed  to  be  specified  functions  of  x,  but  independent  of  t.  In  what  follows  the  gen- 
eral form  of  1 will  be  employed. 

Consider  the  initial  boundary  value  problem  on  xQ  < x^»  for  Equation  47.  At 
x Xq  the  temperature  is  specified 

T (Xq.  0)  = T0  II  (t),  (48) 

where  TQ  is  an  amplitude  constant  and  H (t)  is  the  Heavyside  step  function.  As  x 
it  is  anticipated  on  physical  grounds  that 

I.im  T (x,t)  0 
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with  a similar  condition  on  the  l'lux,  - K HT/  <ix.  The  initial  condition  is  that  no  heat 
is  transferred  prior  to  t 0 or 

T (x,  0)  = 0,  (5°) 

At  this  point  it  is  convenient  to  introduce  dimensionless  variables  as  follows: 

Kn 


11 0 ^O  xo 


X - x/x0,  P = p/pq,  t = t JT 

p, 

T = T/T0,  C C/C0,  K K/Kqi 

where  T0,  C0,  p0,  KQ  are  the  values  of  the  obvious  variables  at  x x0.  Using 
Equations  51,  47-49  and  50  become  respectively. 


(51) 


and 


„ r-  n 

p C 

0 T 

(52) 

<)  X ft  X — 

d t 

T (1,0)  = H(t), 

(53) 

Lim  T (x,t)  0, 

x - 00 

(54) 

T (x,  0)  0. 

(55) 

THE  TRANSFORMED  PROBLEM 

It  will  now  be  shown  that  an  equation  of  the  form  of  Equation  17  is  amenable  to 
analytic  solution  if  the  following  choices  for  K , i>  , and  c arc  made 


K 

O 

X , 

(56) 

p 

= xA, 

(57) 

c 

Xr, 

(58) 

where  u , q and  r are  constants.  Using  Equations  56  to  58,  Equation  is 


( o + 

d x 
0 x _ 


d T “I 

a t 

a + m) 

Y S 

a x — 

A 

a t 

0, 


(59) 
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where 


F 


s = x + r, 

is  introduced  for  convenience. 


(60) 


Thus,  at  this  point,  a solution  to  Equation  59,  subject  to  Equations  53  to  55  is 
sought.  To  this  end  the  Laplace  Transform  of  Equation  59  is  taken,  which,  after  the 
use  of  Equation  55,  is 


' 

d 

dT 

( o + m 

X 

- Xs  p T, 

(61) 

dx  ■ 

dx  — 

In  Equation  61  T is  the  Laplace  transform  of  T and  p is  the  transform  parameter  de- 
fined through 


a(T[(x,t)]  f°°  ePt  T dt  T (x,p), 

J o 


(62) 


In  what  follows,  the  superscript  - will  denote  the  transform  of  the  obvious  quantity. 
The  boundary  conditions,  Equations  53  and  54  take  the  forms,  respectively, 

T (l,p)  = 1/p,  (63) 


Lim  T (x,  p)  0. 
x - 00 


(64) 


It  may  be  recognized,  carrying  out  the  indicated  differentiation,  that  Equation  61 
is  a modified  Bessel  equation,  provided  that  S + 2 / o . The  general  solution  to  that 
equation  is 


T (x,p) 

xa  |j\(p) 

Kn 

(•/) 

■MP>  '„(»  -7). 

where 

a = | 1 

- (m  + a 

) 1 

/2, 

a = Is 

+ 2 - 0 I 

/2, 

n 

a/ 13 

(65) 


(66) 


,1/2 
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The  solution  which  satisfies  the  boundary  condition  63  and  the  initial  conditions, 
but  not  necessarily  the  far  field  condition,  Equation  64,  is 


A(p) 

= [pKn(q/e)]  1 

(67) 

B(p) 

= 0 

(68) 

SHORT  TIME  SOLUTION 


Following  the  method  of  Reference  26  a solution  of  Equations  65,  67  and  68  is 
sought.  Omitting  the  details  (which  are  in  Reference  26),  Equation  65  is  expanded  in 
an  asymptotic  series 


T (x,p)~x 


(a 


-',/2)-i.eXp[--a.(xS-1,]|ltc1/q.e2/q2*...| 


where 


and 


C1  “ A1  " B1 

c2  = A2  ' B2  " B1  (A1  " Bl) 


Al  = -t4"2  :^L  Sx-e 
8 


A2 


(4n2  - 1)  (4n2  - 9)  ^ 2 x~20 
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B, 


(4nz  - 1) 


(69) 

(70) 

(71) 


(72) 


B2  - <4"2  - <4"2  - 9> 

128 


The  inverse  of  the  image  function  given  in  69  may  be  obtained  from  standard 
transform-pair  tables.  The  resulting  approximation  to  T (x,t)  is  then 

T (x,  t)  - x^a  " B//2^-{erfc  (a/Zt1^2)  + Cj  (4t)1//2  i erfc  (a/2t1,/2)  + 

c2  (4t)  i2  erfc  | 


(73) 
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whore 


in  erfc  z 


/°°  n- 1 

i erfc  (,  d $ 

7 


n = 1,2,... 


are  Hartree's  repeated  conplementary  error  functions. 


ACCURATE  APPROXIMATE  SOLUTION 

In  order  to  effect  a more  accurate  solution,  the  method  of  Bickford*29)  is  applied 
to  Equation  65.  In  this  method  the  functions  Kn(z)  are  presented  by 


_ 1 / O 

Kn(z)  - exp  (-z).  [ ” /2Z  J • Pj  (z) 


where 


P\  (z) 


7)  + axZ^  " 1 + . . . aj 

Zi  + bxzj  “ 1 + . . . bj 


is  a rational  function  of  two  polynomials  of  order  j.  The  coefficients  are  chosen  by 
forcing  the  right-hand  side  of  Equation  75  to  yield  exact  values  for  Kn(z)  at  2j  discrete 
points.  For  the  present  paper  adequate  accuracy  was  obtained  by  choosing  j = 1 . The 
values  of  aj,  bj  used  for  this  representation  are  in  Table  A-3.  The  result  of  this  cal- 
culation gives: 

T ft  p1  [(q  + b^  x-0)"1  - (q  + a^)-1]- [p  + q f2  + fal  • exp  [-q  f4]  (77) 


where 


fl  x 


(z  - 1/2)  , , -p.  .. 

' /[a1  - bx  x M ] /3 


f2  = (b4  + ax  x p)P  (78) 

f3  = al  bzp2  X~P 

f4  =j~  <**  - 1) 

The  inversion  of  77  yields 

ft  I/2  l/2 

T (x,t)  ftf5  | exp  [f6V  f6f4]  • erfc  [f6t  + f4/2t  ] 

- exp  [a2o  2t  + aj0f4]  • erfc  [a^t1//2  H f4/2t^2] } + f]^3  erfc[  f4/2t^/2]  (79) 
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TABLE  A-3.  VALUES  OF  at,  bx  USED  IN  BESSEL  FUNCTION  REPRESENTATION 


i 


' 

I 

f 


K i 


> 

i 

i 


al 

bl 

ko 

0.2772 

0.3986 

kl 

0.5696 

0.2023 

kl/4 

0.3031 

0.3949 

kl/2 

1.0 

1.0 

where,  in  addition  to  previously  defined  symbols 
f5  = 0 (b  - a)  (1  - x ,f) 
fg  0bx  ^ 


(80) 


In  general,  the  solution  given  by  Equation  77  is  approximate,  and  the  nature  of 
that  approximation  is  presented  below.  However,  when  k = 0 or  -1  the  solution  ob- 
tained is  exact.  The  instances  of  k = 0 and  -1  correspond  to  homogeneous  ( o = 0)  and 
" 4/3  cases,  respectively.  These  along  with  some  approximate  solutions  are  shown 

in  Figures  A -4  and  A-5. 

Table  A-4  shows  the  values  of  I n I and  p corresponding  to  the  various  values  of 
in  K(x)  = x°  which  are  used  here.  Table  A-3  shows  the  values  of  aj  and  b^  used  in  the 
present  study. 


TABLE  A-4.  VALUES  OF  | n f AND  0 CORRESPONDING  TO  VARIOUS 


o in  k = x°. 

0 

n 

0 

0 

1/2 

1/2 

2/3 

1/4 

2/3 

1 

0 

1/2 

4/3 

1/2 

1/3 

3/2 

1 

1/4 
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Figure  A-4.  Temperature  vs.  Time  at  x = 2 for  Various  Values  of  the  Conductivity:  k(x) 


Figure  A -4  shows  the  ealculation  of  the  temperature,  T as  a function  of  time  at 
x 2 for  n ~ 0,  2/3,  1,  4/3  and  3/2  in  k(x)  = x"  . At  x = 2 the  value  of  the  conductiv- 
ity varies  from  1 (for  <>  0)  to  2.83  (for  o 3/2).  However,  the  temperature  dif- 
ference between  the  « 0 and  <>  - 3/2  curves  is  on  the  order  of  only  25%.  Thus,  even 

with  a very  large  difference  in  conductivity,  the  resulting  temperature  difference  is 
not  great.  It  is  concluded  that,  for  practical  purposes,  one  should  compare  heat  trans- 
fer calculations  based  on  properties  evaluated  at  the  extremes  of  processing  variables. 


IV  FAILURE  MODELS 


It  has  been  tacitly  assumed  in  this  report  that  failure  occurs  only  when  the  mandrel 
punches  through  the  shell,  caused  by  a tensile  failure  in  the  shell  wall.  The  criterion 
for  this  to  occur  was  defined  to  be  the  case  at  which  the  specific  draw  stress  or  draw 
stress  ratio  exceeds  unity.  Thus  it  was  expressed  as 


1 


“"I 

'T  / 


In  this  equation:  is  the  tension  in  the  shell  wall  and  0 is  the  local  yield  stress  in  the 

shell  wall  beyond  the  die  exit.  The  tension,  of,  is  found  as  a result  of  computation 
within  the  model.  It  is  thus  a theoretical  quantity.  The  yield  strength  0 is  determined 
partly  by  theory  and  partly  by  experiment.  As  discussed  earlier  in  this  report,  o is  a 
function  of  temperature,  strain  and  strain  rate.  The  temperature  is  an  assigned 
quantity.  The  strain  and  strain  rates  are  computed  within  the  model.  These  theoret- 
ical quantities  are,  however,  fed  into  experimental  data  when  used  to  compute  the  local 
yield  strength.  The  details  of  this  were  established  on  page  21. 


The  above  failure  model  neglects  the  very  real  possibility  of  plastic  instability. 
Reference  30  discusses  this  concept  in  detail.  Instability  in  drawing  manifests  itself 
by  a large,  sustained,  axial  elongation  or  "necking"  of  the  wall  without  a correspond- 
ing increase  in  applied  drawing  load.  The  point  of  instability,  when  plotted  on  a true 
stress-true  strain  diagram,  is  known  to  be  considerably  lower  than  the  ultimate  (or 
fracture)  strength  of  the  metal,  particularly  in  hot  processing. 

One  important  planned  improvement  to  the  model  is  to  have  the  computer  program 
automatically  test  the  local  yield  strength,  a , for  the  point  of  instability. 


On  the  basis  of  the  above  discussion,  it  is  felt  that  caution  must  be  used  in  design- 
ing dies  and  processing  conditions  which  result  in  a specific  draw  stress  close  to  unity. 
In  such  cases  it  should  be  determined  if  the  instability  limit  has  already  been  exceeded. 
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1.  An  error  has  been  discovered  in  the  theoretical  formulation  and 
computer  solution  for  the  analysis  which  is  given  in  Appendix  A, 
Section  III  on  pages  68  to  76  of  the  subject  report.  As  a result 
the  statement:  "Thus,  even  with  a very  large  difference  in 
conductivity,  the  resulting  temperature  difference  is  not  great", 
which  appears  at  the  top  of  page  76,  is  invalidated.  It  is  further 
advised  that  Appendix  A - Section  III  be  deleted  in  order  to  avoid 
any  possible  misuse. 

2.  The  latter  finding,  however,  in  no  way  invalidates  the  summary 
and  conclusions  of  those  portions  of  the  report  pertaining  directly 
to  artillery  projectiles  (pages  43-46). 


